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Abstract 

We consider the equation 

votB+aB=0 (1) 

on the plane with a being a real-valued function and show that it can 
be reduced to a Vekua equation of a special form. In the case when a 
depends on one Cartesian variable a complete system of exact solutions 
of the Vekua equation and hence of equation (1) is constructed based 
on L. Bers' theory of pseudoanalytic formal powers. 



1 Introduction 

Solutions of the equation 

rot B + a B = (2) 
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where a is a scalar function of space coordinates are known as Beltrami 
fields and are of fundamental importance in different branches of modern 
physics (see, e.g., [22], [18], [7], [2T], [I], [9], [8], [II]). For simplicity, in this 
work we consider the real-valued proportionality factor a and real-valued 
solutions of (T5]), though the presented approach is applicable in a complex- 
valued situation as well with a considerable complication of mathematical 
techniques involved (instead of complex Vekua equations their bicomplex 
generalizations should be considered [5] , [H] ) • We consider equation Q ona 
plane of the variables x and y, that is a and B are functions of two Cartesian 
variables only. In this case as we show in section [3] equation (T5]) reduces to 
the equation 

div Q Vm ) +au = 0. (3) 

This second-order equation can be reduced (see [H]) to a corresponding 
Vekua equation (describing generalized analytic functions) of a special form. 
This reduction under quite general conditions allows us to construct a com- 
plete system of exact solutions of ([3]) explicitly (see [13] and [IS]). For the 
reduction of ([3]) to a Vekua equation it is sufficient to find a particular solu- 
tion of ([3]). In the present work (section H]) we show that in a very important 
for applications case of a being a function of one Cartesian variable a particu- 
lar solution of ([3]) is always available in a simple explicit form. This situation 
corresponds to models describing waves propagating in stratified media (see, 
e.g., [IS])- As a result in this case we are able to construct a complete system 
of solutions explicitly which for many purposes means a general solution. We 
give an example of such construction. 

We show in this work that when a = a(y) (of course in a similar way 
the case a = a(x) can be considered) equation ([3D and hence equation ([2]) 
reduce to the Vekua equation of the following form 

d*W{z,y) = t ^W(z,y) (4) 

where 

/ = — = sin A H — = cos A; 

'a \ a 



A is an antiderivative of a with respect to y, c\ and C2 are arbitrary real con- 
stants, z = x + iy and = \{d x + id y ). A complete (in a compact uniform 
convergence topology) system of exact solutions to (j4j) can be constructed 
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explicitly. The system represents a set of formal powers [3J, [B] which gen- 
eralize the usual analytic complex powers (z — z ) n , n = 0, 1, 2, . . . and in 
a sense give us a general solution of (j4j). Thus, in the case when a is a 
function of one Cartesian variable the Vekua equation equivalent to (j2J) in a 
two-dimensional situation can be solved and a complete system of solutions 
of fl2j) is obtained. 



2 Preliminaries 

We need the following definition. Consider the equation 

d& = $ (5) 

on a whole complex plane or on a convex domain, where $ = $i + i$2 is a 
given complex valued function such that its real part $i and imaginary part 
$2 satisfy the equation 

dy^ - d x § 2 = (6) 

then there exist real valued solutions of (jHJ) which can be easily constructed 
in the following way 

ip(x,y) = 2U ^ 1 (rj,y)drj + /% 2 (x ,£R) +c (7) 

where (xo,yo) is an arbitrary fixed point in the domain of interest and c is 
an arbitrary real constant. 

By A we denote the integral operator in ([7]): 

A[$\(x,y) = 2[( 1> 1 ( V ,y)d v + f" $ 2 (x , £)(%)+ c. 

\Jxo Jyo J 

Note that formula ((7|) can be extended to any simply connected domain by 
considering the integral along an arbitrary rectifiable curve V leading from 

(x , y ) to (a;, y) 



(p{x, y) = 2[ / §xdx + § 2 dy + c. 



r 



Thus if $ satisfies (Q, there exists a family of real valued functions <p such 
that dzf = $, given by the formula ip = A[<&]. 
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Let / denote a given positive twice continuously differentiable function 
denned on a domain Q G C Consider the following Vekua equation 



= jW in n (8) 

where the subindex z means the application of the operator 8%, W is a 
complex- valued function and W is its complex conjugate function. As was 
shown in [12] , [13] , [H] , [15] , [TF] equation (jSJ) is closely related to the second- 
order equation of the form 

(divpgrad +q)u = in Q (9) 

where p and q are real-valued functions. In particular the following state- 
ments are valid. 

Theorem 1 [L$ Let u be a positive solution of (TJ|). Assume that f = ^/pu 
and W is any solution of (TJj]. Then u = ^=ReW is a solution of (TJ|) and 
v = y/p Im W is a solution of 

(div - grad +qi)v = in Q (10) 
V 



where 



, 1 = _i(i +2 ^,^ +2 f^V|. 



V \ V \ P u I \ u 



Theorem 2 [Lffl Let Q be a simply connected domain, Uq be a positive so- 
lution of (Tj|) and f = p 1 / 2 ti . Assume that u is a solution of (Tj|). Then a 
solution v of pop with q\ defined by ( TO)) such that W = p l l 2 u + ip^^v is a 
solution of is constructed according to the formula 

v = u^A^pul^u^u)) (12) 

and vice versa, let v be a solution of (fiOJ). then the corresponding solution 
u of (TjJ) such that W = p x l 2 u + ip~ l l 2 v is a solution of is constructed 
according to the formula 

u = — u$A{ip~ l UQ 2 dz(uov )). (13) 
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Thus the relation between (jS]) and is very similar to that between 
the Cauchy-Riemann system and the Laplace equation. Moreover, choosing 
p = 1, q = and «o = 1 we obtain that ffT2l) and f|T3|) become the well known 
formulas from the classical complex analysis for constructing conjugate har- 
monic functions. 

For a Vekua equation of the form 

Wz — aW + bW 

where a and b are arbitrary complex-valued functions from an appropriate 
function space [20J a well developed theory of Taylor and Laurent series in 
formal powers was created (see [3J, [I]) containing among others the expan- 
sion and the Runge theorems as well as more precise convergence results 
(see, e.g., [19]) and a general simple algorithm [15] for explicit construction 
of formal powers for the Vekua equation of the form ([8]). 



3 Reduction of (GS) to a Vekua equation 

We consider equation ([2]) where both a and B are supposed to be dependent 
on two Cartesian variables x and y. Then equation ([2]) can be written as the 
following system 

d y B 3 + aB 1= (14) 
- d x B 3 + aB 2 = (15) 
d x B 2 - d y B 1 + aB 3 = 0. 
Solving this system for B% leads to the equation 

A5 3 - VB 3 ^ + a 2 B 3 = (16) 

where (■, ■) denotes the usual scalar product of two vectors. 
Note that 

adiv ( -VS 3 ) = AB 3 - ( — , V5 3 

\a J \ a 

and hence (1T61) can be rewritten as follows 

div (Iv5 3 )+c*53 = 0. (17) 
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Thus equation fl2]) reduces to an equation of the form (P) with p = 
q = a. 

Let us notice that (see, e.g., [14"] ) 



\ja and 



div -V + a = —= (A - r) -= 

where 

IAq 3 fVa\ 2 2 , . 

r = + - -a 2 . 18 

2 a 4 V a / 

That is B 3 is a solution of f[TTj) iff the function / = B 3 /y/a is a solution of 
the stationary Schrodinger equation 

(-A + r)/ = (19) 

with r defined by ([IB"]) . As was explained in section [21 given its particular 
solution this equation reduces to the Vekua equation (jSJ). Unfortunately, in 
general we are not able to propose a particular solution of ([17]) . Nevertheless 
in an important special case when a depends on one Cartesian variable, a 
particular solution of ([IT]) is always available in explicit form. We give this 
result in the next section. 



4 Solution in the case when a is a function of 
one Cartesian variable 

Let us consider equation f[T9]) where a = a{y). We assume that a is a 
nonvanishing function and look for a solution of the corresponding ordinary 
differential equation 



dy 2 12 a 4 \ a 

Its general solution is known (see [101 2.162 (14)]) and is given by the ex- 
pression 

f (y) = -p=sm A{y) + -^=cos A(y) (20) 
where A is an antiderivative of a and cj, C2 are arbitrary real constants. 
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Choosing, e.g., c\ = 1, c 2 = and calculating the coefficient (dzfo) / fo 
we arrive at the following Vekua equation which is equivalent to (j2J) in the 
case under consideration (and which is considered in any simply connected 
domain where sin A(y) does not vanish), 

dzW(x,y) = % - (a(y)zoXA{y) - W(x,y). 

Note that F = f = sm -^} and G = 4- = ^"ffl represent a generating 

pair for this Vekua equation (see [13], [15]) and hence if W is its solution, 
the corresponding pseudoanalytic function of the second kind u = i Re W + 
ifo Im W satisfies the equation 



= T-rifis ( 21 ) 



l + /o 2 

which can be written in the form of the following system 



rtYyi 4>y r-2 

Jo Jo 



where = Re a; and ip — Imw. 

For /o being representable in a separable form f (x,y) = X(x)Y(y) the 
formulas for constructing corresponding formal powers explicitly were pre- 
sented already by L. Bers and A. Gelbart (see [3] and [6]). Using them we 
obtain the following representation for the formal powers corresponding to 

m 

*Z^(a,z ;z) = ai J2 (l) {x-x f n ~ k \ k Y k 

k=0 ^ ' 

+ ia 2 ^( n k \{x-x^ n ~ k \ k Y k 

k=0 ^ ' 



(we preserve the notations from [3J) where z Q = x^ + iyo is an arbitrary point 
of the domain of interest, a is an arbitrary complex number: a = d\ + ia 2 , 
Y k and Y k are constructed as follows 

r(°W) = ? (0 W) = i 
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and for n = 1, 2, . . . 

KW(y Q ,y) = n I Y^iyo^fS^dv n odd 



y 

Y^(y ,y)=n[Y^ 1 \y ,r ] ) 



n even 



yo 

y 



Y^(x ,x) = n lY^iyo^)-^- n odd 



yo 
v 



Yto(x ,x) = n / Y^Xy^fMdr) n 



even. 



yo 



The system {*Z( n \l, zo', z), *Z( n \i, z ; z)} _ represents a complete (in a 
compact uniform convergence topology [2]) system of solutions of (}2~Tj) that 
means that any solution u of (1211) in a simply connected domain Q can be 
represented as a series 



oo 

UJ{Z) = 

n=0 n=0 



£ a {n \a n , z ; z) = J2 K * Z(n) (l, ^o; z) + < ^ (n) (^ ^o; z)) 



where a' n = Rea n , = Ima„ and the series converges normally (uniformly 
on any compact subset of fl). Consequently the system of functions 

{/o(y) Re(*Z^(l, z Q ; «)), /o(y) Re(*Z (r % z ; ^))}7 =0 

represents in the same sense a complete system of solutions of (1191) with r 
defined by (|T8|) . and 

{ \ffiy)Mv) Re(^ (n) (i, z ; *)), Vrty)fo(y) M*z {n) (i, z 0] z))}°° (22) 

l. J n=0 

is a complete system of solutions of (|17[) . Thus in the case under consideration 
any solution B 3 of (|17|) can be represented in the form 

00 

B 3 (x, y) = ^(a n sin A{y) Re{^ n \l, z ; z)) + b n wxA{y) Ee(^ (n) (z, z ; z))) 

n=0 



where a n and b n are real constants. 

The other two components of the vector B are obtained from (1141) and 
(US]): 

B x = d v B 3 and B 2 = -d x B 3 (23) 

a a 

that gives us a complete system of solutions of ([2]) in the case under consid- 
eration. On the following example we explain how this procedure works. 

Example 3 Let us consider the following relatively simple situation in which 
the corresponding integrals are not difficult to evaluate. Let 

«(y) = -t== 2 ( 24 ) 

and Q be an open unitary disk with a center in the origin. We take in l[2U\) 
c\ — and c 2 = 1 . Then it is easy to verify that 

f (y) = (l-y*)i. 

The first three formal powers with a centre in the origin can be calculated as 
follows 

x r ?/(l - ?/ 2 )^ 3y(l-y 2 )^ 3 

*Z«(1, 0; z) = x + ? P A y ' + \ + g arcsiny], 



*zW(i,0;z) = - V -3 +ix, 

(i -y 2 ) 2 

Z^il 0-z) = x 2 - 1 2_3 yaxcsiny 
* 1 ,U ' j X 4^ 4(i_ y 2)i 



, y(l-y 2 )l 3y(l-y 2 )^ 3 . 
+ 22X I 1 h - arcsm y ] , 



^ (2) (<, 0; *) = —^-r + i(x 2 -y 2 -\y* 



(l-^f 



2* 
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^>(l,0;,)=x 3 -3x( l y 2 + 3 yaxcsiny 
1 ^ U 4(i_j,2)i 



2/ ^(1 -y 2 )§ , 3y(l -y 2 )^ , 3 . \ 
+ 3zx I 1 h - arcsm y I 

- <(-^y(i - y 2 ) 1 + ^y(i - y 2 ) 1 + y(i - y 2 )^ - ^y 2 ) 

— — fl — y 2 ) 2 arcsiny + — — arcsiny), 
lg v y ) v 12g y;> 

^ 3 )(*,0;z) = 3 y(l + y 2 ) 3 / 2 _ / 2 _ 1 A\ 

l ' J 4(1-^)1 4 v v 2 y 

Now taking the real parts of these formal powers and multiplying them by the 
factor y/afo (see §EBj) ) we obtain the first elements of the complete system 
of solutions of (T7\ ), that is any solution B 3 of fcFfy in a simply connected 
domain can be represented as an infinite linear combination of the functions 

2 i ( 2 . 2-/2 1 2 ^ y arcsin y 



{(1-y 2 )^, x(l-y 2 )5, -y, (1 — 2/^)2 \x 2 --y 



4 4(1-3,2)1 



-2 Xy , i i- yr , { ^-3 X[l f+-y-—-y 

- 3x 2 y + -y(l + y 2 ) - -(1 - y 2 )^ arcsiny, ...} 

and the corresponding series converges normally. 

From ( d^ ) it zs easy to calculate the corresponding components B\ and B2 
respectively, 

3 3 5 

{y, xy, (l-y 2 )^, - (1 -y 2 ) I arcsin y + y(:r 2 - -y 2 + -), 

2x(l — y 2 ) 5 , ~~ y 2 ) 5 arcsiny + y ^a: 3 — ~^y 2 + — x 

- (^y 2 - 3x 2 )(l -y 2 )5 - ^y arcsiny, ...} 

and 

{0, (1-y 2 ), 0, 2a;(l-y 2 ), -2y(l-y 2 )^ 
(3x 2 - ^y 2 )(l - y 2 ) - ^y(l - y 2 )^ arcsiny, -6xy(l - y 2 )^, ...}. 
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Thus, we obtain the following complete system of solutions of (TJ|) with the 
proportionality factor a defined by 



B, 



y 
o 



y 2 ) 2 



xy 

(i-y 2 ) i 

x(l — y 2 )2 



B. 



1 if) 



-y 



B- 



I 2 ) 2 arcsin y + y(x 2 — \y 2 + f ) \ 



2x{\ - y 2 



y 2 ) 2 



„2 1„,2 3 y arcsin y 
x ~ jy ~ I~ 7~T 



B 



2x{l-y 2 )^ 

-M^-y 2 ) 

—2xy 



B, 



/ — y 2 ) 2 arcsin y + y (x 3 — \xy 2 + ^fx) \ 
(3x 2 - h 2 ){l - y 2 ) — f?/(l — y 2 )^ arcsin y 



\ 



y 2 Y 2 



:h 2 



x 



Q T / 1„2 i 3 y arcsin y \ \ 



3a; 2 ) (1 — y 2 )z — jy arcsin y 



B, 



—3x 2 y 



-6xy(l — y ) 2 
H 



h(i + y 2 )-Hi-y 2 ) 



arcsm y 



5 Concluding remarks 

We presented two new results. 

1. We showed that the system of equations describing Beltrami fields on 
the plane can be reduced to a Vekua equation of the form §8§ whenever 
any particular solution of the corresponding second-order equation of 
the form (fl7|) or (fl9|) is known. 

2. In the case when the proportionality factor a depends on one Carte- 
sian variable we obtain a particular solution explicitly, construct the 
corresponding Vekua equation (section 4), and solve it in the sense that 
a complete system of solutions is obtained which gives us a complete 
system of solutions of the original vector equation describing Beltrami 
fields. 
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Of course not always the integrals involved in the construction of the 
complete system of solutions are sufficiently easy to evaluate explicitly as 
in the example [3j Nevertheless our numerical experiments confirm that in 
general the formal powers and hence the solutions of (j2J) can be calculated 
with a remarkable accuracy. For example, the vector B ^ (see notations in 
the example [3]) in the Matlab 7 package on a usual PC can be calculated with 
a precision of the order 10~ 4 . Thus, the use of formal powers for numerical 
solution of boundary value problems corresponding to (121) and more generally 
to equations of the form ([9]) is really promising. The work in this direction 
will be reported elsewhere. 
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